from scipy.linalg import solve
A = np.array([[1, 2], [3, 4]])
b = np.array([5, 6])
solution = solve(A, b)
print(f"Solving linear system: {solution}")